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Various models of QCD vacuum predict that it is dominated by excitations that are predominantly 
self-dual or anti-self-dual. In this work we look at the tendency for self-duality in the case of 
pure-glue SU(3) gauge theory using the overlap-based definition of the field-strength tensor. To 
gauge this property, we use the absolute X-distribution method which is designed to quantify the 
dynamical tendency for polarization for arbitrary random variables that can be decomposed in a 
pair of orthogonal subspaces. 
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1. Motivation 

Various models of QCD vacuum use semi-classical arguments to describe the mechanism re- 
sponsible for confinement or chiral- symmetry breaking. The semi-classical arguments start by 
expanding QCD partition function around extremal points of the action, i.e., 
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The first task is then to find the extremal points of the action and then take into account gaussian 
fluctuations around these extrema. 

The action for pure-glue QCD can be expressed in terms of the self-dual and anti-self-dual 
components of the field strength tensor 
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The integral of the F^^F^y term is a boundary term that is related to the topological charge of the 
configuration. If we keep the boundary values fixed, the integral is minimized when the quantity 
in the parenthesis vanishes. This happens when the field is self-dual, F^^ = F^y, or anti-self-dual 
F^^, = —F^y. A more sophisticated analysis leads to the conclusion that all the extremal points of 
the classical action that are not saddle points satisfy this condition [1]. 

It is then natural to expect that if QCD vacuum is correctly described by a semi-classical model, 
the field strength in a typical lattice QCD ensemble will exhibit a high degree of self-duality. To 
gauge this tendency we decompose the field strength at every point on the lattice into its self-dual 
components and analyze their polarization properties. To do this, we use the method of absolute X- 
distribution designed to analyze the dynamical aspects of polarization [2]. A more detailed account 
of this work is given in Ref. [3]. 



2. Dynamical polarization 

We start by reviewing the method of absolute X-distribution. A first version of this approach 
was introduced in a study of the local chirality of the low-lying eigenmodes of the Dirac opera- 
tor [4]. In general, for an arbitrary observable that can be split in two components Q = Q\+ Qi, we 
say that Q is polarized when it tends to be aligned with either one of the components. More pre- 
cisely, if we look at the magnitude of components, qi= \\Qi\\, we tend to think that the observable 
Q is polarized when the probability distribution J^b{Qi,(l2), with support in the positive quadrant 
of the qiq2-plane, is peaked in the vicinity of the <7i 2 axes. 

The raw distribution ^b{<li,<l2) is difficult to characterize. A more direct measure is offered 
by the induced distribution of the polarization angle. In Fig. 1 we plot the raw distribution of 
chirality components as determined in a previous study [2] and the corresponding polarization 
angle distribution (the curve indicated by a = 1), which we call the X-distribution. We see that 
the X-distribution tends to be concentrated towards the middle of the graph, suggesting an anti- 
polarization tendency. 
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Figure 1: Sample pair distribution generated by chirality components of the lowest eigenmodes of ensemble 
£1 from [2]. Right: the associated X-distribution, i.e., the induced distribution of the polarization angle. 

A more careful analysis reveals that the conclusions based on this method can be misleading. 
The X-distribution is determined by the choice of parametrization for the angles measured in the 
^1^2-plane. The definition we used to plot Fig. 1 is 

4 IIG2II 

x = — arctan rr — 1 . (2. 1) 
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We will refer to this choice as the reference polarization [4]. However, this choice is not unique. 
Alternative definitions were used in various studies. Using t = ||22||/||2i ||, one class of valid angle 
variables is given by a generalization of the above definition 

4 

f = -arctan(f")-l, (2.2) 
% 

where a > is an arbitrary parameter [2]. For a = 1 the angle parameter x is the reference 
polarization defined above, while the definition based on x with a = 2 was used in a study of self- 
duality in pure gauge QCD [5]. In the right panel of Fig. 1 we compare the X-distribution for the 
ensemble shown in the left panel, measured using the reference polarization and the polarization 
defined by x with a = 4. The qualitative behavior of the distribution changes dramatically, while 
the dynamics producing the original distribution is unchanged. It is clear then that conclusions 
based on X-distributions alone cannot be trusted. 

To address this problem we define the absolute X-distribution, a measure of the pair correlation 
induced by the underlying dynamics [2, 6]. The basic idea is to compare the correlated distribu- 
tion ^b{(]i,Q2) with a similar distribution where the components are statistically independent, to 
isolate the effect of the dynamics. The uncoiTclated distribution is constructed from the marginal 
distributions 

P\{q\) = j dq2 3^b{qi,q2) and Piiqi) = j dqi 0^h{q\,qi) ■ (2.3) 

For our application, symmetry guarantees that Pi = P2- The uncorrected distribution is S^u{q\ , qi) = 
P\{q\)P2{q2)- We define an angle variable that has constant angular density for the uncorrelated 
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Figure 2: X-distribution using the reference polarization for the correlated and uncorrelated distributions 
(left) and the absolute X-distribution (right). 



distribution. This is the absolute polarization. The histogram of this angle variable for the un- 
correlated distribution is flat. In our figures this is indicated by a horizontal dashed line. The 
X-distribution in terms of the absolute polarization is the absolute X-distribution. 

In the left panel of Fig. 2 we present the X-distribution for the reference polarizations for both 
correlated distribution, I^b, and the uncorrelated one, for the ensemble presented in Fig. 1. 
Notice that these two distributions are almost identical indicating that there is little dynamical 
correlation. In the right panel we plot the absolute polarization histogram, which is almost flat. 
There is a small enhancement towards the edges indicating that the dynamics induces a slight 
polarization. This is consistent with the plots in the right panel, where we see that the uncorrelated 
distribution is more prominent towards the center of the histogram. 

Based on the absolute polarization distribution, Pa{x), we construct a more compact measure 
of the polarization tendency, the correlation coefficient 

1 

Ca = 2r-l where T = j dxPA{x)\x\ . (2.4) 

-1 

The coefficient Y measures the probability that a sample drawn from distribution is more po- 
larized than one drawn from When we have no dynamical correlation this probability is 0.5; 
the correlation coefficient is scaled such that Ca = in this case. 

3. Field strength definition 

In this study, we will use a definition of the field strength based on the overlap operator. 
Compared to the ultra-local definitions, the overlap definition is less susceptible to ultra-violet 
fluctuations, so no arbitrary link smearing or cooling is needed. Moreover, this definition provides 
a natural expansion in terms of eigenmodes of the Dirac operator which allows us to define a 
smoothed version of field strength tensor controlled by the value of the eigenvalue cutoff. 

If we denote with Sp = V/D(.x,3;) the fermionic contribution to the action in the overlap for- 
mulation, it is easy to show that tr^ OfivD{x,x) has the same quantum numbers as the field strength 
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Ffxv [7]. Here ti\ denotes the trace over the spinor index. It was shown by expUcit calculation that 
on smooth fields in the limit a — )• these definitions agree [8, 9], 

tr,a^v£>(x,x) = JF^y{x) + ff{a^) . (3.1) 

Above, is a constant that depends on the kernel used to define the overlap operator. The lattice 
version of the field strength operator used in this study is 

-^MvW = -^tr,o^vD{x,x) = -^tVs0^v[2p-D{x,x)] , (3.2) 

where 2p is the largest eigenvalue of D, the eigenvalue associated with the zero modes' partners. 
We used the fact that tr^ a^v = to cast the definition in a form useful for eigenmode expansion. 

Using the expansion in terms of the eigenmodes of the Dirac operator, we define the smoothed 
version of the field strength [2] 

fI^^{x) = -\ L tr,a^v(2p-A)i//AWV^AW^ (3-3) 

|A|<Aa 

This definition has the property that limA^«. = F and that the contribution of the largest eigen- 
modes is suppressed. The self-dual and anti-self-dual parts of the field strength are defined using 
the dual of the field strength y = ^^iivapFap 

Fs = \{F+F) Fa = \{F-F). (3.4) 

4. Numerical results 

For our study we used a set of pure-glue ensembles generated using Iwasaki action [10]. The 
parameters for these ensembles are presented in Table 1 . To study the continuum limit we have a 
set of 5 ensembles with the same volume. To determine the finite volume effects we also generated 
one ensemble with a larger volume. 

In Fig. 3 we plot the histogram for the absolute polarization for all ensembles with volume 
(1.32fm)'*. We find a small tendency for polarization that decreases as we make the lattice spacing 
smaller. To understand whether this tendency survives the continuum limit, we compute the cor- 
relation coefficient and fit it with a quadratic polynomial in a. As we can see from the right panel 
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Table 1: The size and lattice spacing for the ensembles used in this study. 
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Figure 3: Left: absolute X-distribution for self-duality components. Note that the y-scale is magnified to 
better show the difference between different lattice spacings. Right: the correlation coefficient as a function 
of the lattice spacing and its continuum hmit extrapolation. Error bars are present in these plots but they are 
smaller than the symbol size. 



of Fig. 3 the polynomial fits the data well. The coefficient remains positive in the continuum limit, 
indicating a very small tendency for polarization. The probability that the sample drawn from the 
correlated distribution is more polarized than one drawn from the uncorrelated distribution is 5 1 % 
compared to 50% when the dynamics would produce no correlation. 

To gauge the size of the finite volume effects, we compute the absolute polarization on two 
ensembles with the same lattice spacings but different volumes. Referring to Table 1, these are 
ensembles £"4 and E^. In the left panel of Fig. 4 we compare the absolute polarizations on these 
two ensembles. We find no difference between the two histograms and we conclude that the finite 
volume effects are negligible. 

We also computed a set of eigenmodes of the overlap Dirac operators on ensembles E2, £3 
and £4 and used them to compute the smoothed field strength operator F^. To study the continuum 
limit, a consistent definition of the smoothed operator sums over all modes smaller than a physical 
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Figure 4: Left: absolute X-distribution for ensemble E4 (circles) and E(, (crosses) which have the same 
lattice spacing but different volume. Right: correlation coefficient for the smoothed strength field (diamonds) 
compared to the full version (circles). Error bars are included in both plots. 
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Figure 5: Left: X-distribution for self-duality components of a smooth field strength based on the low-lying 
modes of the chirally-improved Dirac operator [5]. The curved marked with 100% is the relevant one for 
our comparison. Right: absolute X-distribution Pa and X-distribution P, based on two different polarization 
variables (see Eq. 2.2) for ensemble El- 



cutoff. We set the cutoff A = lOOOMeV and found that the behavior of the absolute X-distribution 
is similar to the full version of the operator. In the right panel of Fig. 4 we compare the correlation 
coefficient with the one computed using the full operator. We find that while the values of the 
correlation coefficient ai^e slightly different, the qualitative behavior remains the same. 

We conclude our discussion with a comparison with a similar work by Gattringer [5] who 
studied the self-duality polarization using a smoothed field strength operator. This operator was 
constructed using an eigenmode expansion of the chirally-improved Dirac operator. In Ref. [5] it 
was found that the self-duality exhibits a strong polarization (see left panel of Fig. 5) supporting a 
model of vacuum dominated by topological "lumps". In contrast, we only find a mild dynamical 
tendency for polarization. This is seen in the right panel of Fig. 5 where we plot the absolute X- 
distribution of ensemble E2 which is similar to the ensemble used in Ref. [5]. The discrepancy 
is due to the fact that Ref. [5] uses a polarization measure dominated by kinematical effects. To 
show this, in the right panel of Fig. 5 we also plot the X-distribution measured using the reference 
polarization, a = 1, and the polarization angle used in Ref. [5], a = 2. To better compare our 
results, for these plots we used, as in the referenced study, a smoothed constructed using the 
same number of modes. We see then that when using the same angle definition, our results are 
consistent with those of Ref. [5]. However, using another valid angle parametrization produces 
qualitatively different results due to kinematical effects. We conclude that the strong polarization 
observed in Ref. [5] is mainly due to the specific choice of angle variable rather than the underlying 
dynamics. 

5. Conclusions 

In this work we studied the dynamical polarization propertied of self-duality components in- 
duced by pure-glue QCD dynamics. We found a very mild polarization tendency that survives in the 
continuum limit. This results has negligible finite-volume corrections. The self-duality tendency 
is very small making it unlikely that the vacuum fluctuations are well-described by semi-classical 
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models. Our findings are at variance with the results of a previous study [5]. We conclude that the 
discrepancy is the result of kinematical effects. 
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